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Abstract 

J> ' We consider a uniform distribution on the set A4fc of moments of order k G 'N corresponding 

^ ' to probability measures on the interval [0, 1]. To each (random) vector of moments in 7W2n-i 

we consider the corresponding uniquely determined monic (random) orthogonal polynomial 
^ ' of degree n and study the asymptotic properties of its roots if n ^ oo. 

o 

g ; AMS Subject Classification: 60F15, 33C45, 44A60 

\ Keyvv^ord and phrases: moment space, random moment sequence, random orthogonal polynomial, 
■ ^'^c^i'^^ distribution, Chebyshev polynomials, random matrices. 

1 Introduction 

For a probability measure r] on the interval [0, 1] let 

(1.1) Ck = Ck{v)= [ x''r]{dx) k = 0,1,2,... 

Jo 

denote the corresponding moments and consider the moment space (of order n) defined by 

(1.2) = {{cii^r]) , . . . , Cn{i]))^ I ?7 is a probability measure on the interval [0, 1] } . 

2 

The set Ain is a very small subset of R" with volume proportional to 2~" and has been studied 
extensively in the literature [see e.g. Karlin and Shapeley (1953), Skibinsky (1967, 1968, 1969), 
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Chang, Kemperman and Studden (1993) or Gamboa and Lozada-Chang (2004)]. In order to 
understand the structure of the moment space Chang, Kemperman and Studden (1993) assigned 
a uniform distribution over and studied the asymptotic properties of the first k components 
of the random vector 

(1.3) {Ci,n,...,Cn,nf r.U{Mn) 

if n — oo. In particular, they showed that an appropriately standardized version of the vector 
(Ci,„ . . . , Ck,nY converges weakly to a centered normal distribution, that is 

(1.4) {(Ci,„, . . . , C,,„)^ - (c?, . . . , clY) ^ A^(0, E), 
where denotes the kih moment of the arcsine distribution defined by 

and the matrix E e M*^^*^ is given by 

(1-6) S=i(e°«-e»c°):-,^^, . 

A large deviation principle for the random moment vector (Ci_„ . . . , Ck^rii^ was derived in Gamboa 
and Lozada-Chang (2004), while Dette and Gamboa (2007) investigated the asymptotic properties 
of a moment range process. 

It is the purpose of the present paper to provide further insight in the probabilistic proper- 
ties of quantities associated to random moment sequences. In particular, by the definition 
C{x^) = Ck,2n-i (/c = 1, . . . , 2n - 1) each random vector (Ci,2n-i, ■ ■ ■ , C2n-\,2n-iY -^U{M2n-i) 
defines a (random) moment functional and a corresponding sequence of monic random orthogonal 
polynomials Po,n{x), Pn,n{x) satisfying 

(1.7) £(Pfe,„ P;,„) = if k^l - k,le{Q,...,n} 

[see Chihara (1978)]. These polynomials are carefully introduced in Section 2, where we also 
state some non-standard results regarding moment theory. In Section 3 we derive the asymptotic 
properties of the (random) roots of the polynomial Pn,n{x) associated with the 

random moment sequence (Ci,2n-i5 . . . , C'2n-i,2n-i) ~ W(-^2n-i). In particular, it is shown that 
the empirical distribution function of the (random) roots . . . „ converges almost surely 
to the arcsine distribution if n — cxd and that an appropriate standardization of the vector 
Xi^ni ■ ■ ■ 7 Xk^n is asymptotically normal distributed, where the roots of the Chebyshev polynomial 
of the first kind are used for the centering, and the rate of convergence is 1/ \/n. 
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2 Random coefficients in a three term recurrence relation 



The set A4n defined by p.2p is a very small compact subset of the unit cube [0, 1]" with nonempty 
interior and volume 



(2.1) 



[see Karlin and Shapely (1953)]. The interior of Ain is denoted by throughout this paper. 
It is well known that there exist an infinite number of probability measures on the interval [0, 1] 
with moments up to the order 2n — 1 given by (ci, . . . , C2n-i _i [see Dette and Studden 

(1997)]. Moreover, there exists a unique measure ri~ supported on exactly n points in the open 
interval (0,1) such that Ci{rj~) = J^x^dri~ (x) = Ci{i = l,...,2n — 1) and such that the point 
(ci(r/~), . . . , C2n{v~))'^ ^ dAi2n whcrc 9A^2n dcuotcs the boundary of the set A^2n- The measure 
r]~ is called lower principal representation of the point (ci, . . . , C2n-i)^ € A^2n-i l^^e Skibinsky 
(1967)]. A straightforward calculation shows [see e.g. Szego (1975) or Chihara (1978)] that for a 
given vector 

(Ci, . . . C2„-l)^ = (Ci(?7), . . . , C2n-liv)f e A<2n-1 

the polynomials Po,n(a^) = 1, 



(2.2) 



Co ■ ■ ■ Cm-1 1 



C2m-1 X, ' 



Co 



Cm— 1 



Cm— 1 



C2m-2 



m 



n 



have leading coefficient 1 and are orthogonal with respect to the measure rj , that is 



(2.3) 



Pk,n{^)PiA^)dv (x) =0 if k I, 



where the moment functional £ is defined by 

• 1 



(2.4) 



C(x'')=Ck= I x^d'q-^x) (A; = 0,...,2n- 1) . 





In other words: each vector (ci, . . . , C2„--i)^ G A^2n-i uniquely determines monic orthogonal poly- 
nomials Po,n(a;), . . . ,P„,n(a;) satisfying (12. 3p . In the following we consider a one-to-one mapping 
of the set M.2n-i on the cube (0,1)^""^, which was introduced by Skibinsky (1967). The new 
coordinates are called canonical moments and have been studied by numerous authors [see e.g. 
Dette and Studden (1997) for a detailed discussion]. To be precise, let ^([0, 1]) denote the set of 



3 



all probability measures on the interval [0, 1], ^k{x) = (x, . . . , x^) the vector of all monomials of 
order k and define for a fixed vector c = (ci, . . . , Ck)^ G J^k 



Sk{c) := |/i G V{[0, 1]) : ^ <l>k{x)fi{dx) = c| 



as the set of all probability measures on the interval [0, 1] whose moments up to the order k 
coincide with c = (ci, . . . , Cfc)^. Note that this set is a singleton if and only if c G dJ^k [see Dette 
and Studden (1997)]. For k = 2,3, .. . and for a given point (ci, . . . ,Cfc-i)^ G A4k~i we define 
= c^(ci, . . . , Cfc_i) and = c^(ci, . . . , Ck-i) as the largest and smallest value of Ck such that 



ci, 



. , Ck)^ G 9A^fc, that is 



= mm 



max 



x^lJL{dx) I /i G S'fc-i(ci, . . . , Cfc-i)|, 
x^^{dx) I ;U G 5'fc__i(ci, . . . ,Cfc_i)|. 



Note that < Ck < and that both inequalities are strict if and only if (ci, . . . , Ck-i) G A^^^i 
[see Dette and Studden (1997)]. For a moment point c = (ci, . . . , c„)^ in the interior of the moment 
space A^„ the canonical moments or canonical coordinates of the vector c are defined by 

(2.5) pi = ci, and pk = — ^ k = 2, . . . ,n . 

Ck -Ck 

Note that < < 1 (A; = 1, . . . , tt.) if (ci, . . . , c„)^ G , and that the definition (12. Sp defines 
a one to one mapping between and the open unit cube (0, 1)". For more details regarding 
canonical moments we refer to the work of Skibinsky (1967,1968,1969) and to the monograph of 
Dette and Studden (1997). In particular it is shown in the lastnamed reference that the three 
term recurrence relation corresponding to the monic orthogonal polynomials defined by (12. 2p can 
be represented in terms of canonical moments, that is Po{x) = 1; Pi{x) = x — (i , and for 
1 < m < n — 1 

(2.6) Pm+l{x) = {x — ~ (2m+l)Pm{x) — (2m-lC2mPm-l{x), 

where the quantities Cfe are given by Co = 0, Ci = Pii Cfc = Ik-iVk with g^-i = 1 — Vk-i if ^ > 2. 
Moreover, from the representation ( 12. 5p it is easy to see that dci/dpi = 1, 

'^Pj [ rk{ci, . . . ,Cfc_i) if J = k 

{k > 1) where 



k 



(2.8) rk+i{ci, . . . , Cfc) = Tfc+i = - C;^^^ = Y[pj{l - Pj) 
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denotes the range of the moment space A4k (with the convention ri = 1). In the following sec- 
tion we will use these results to study stochastic properties of the roots of random orthogonal 
polynomials associated with a uniform distribution on the moment space A^2n-i- 

3 Asymptotic zero distribution of random orthogonal poly- 
nomials 

For each G N let C2n-i = (C'i,2n-i5 • • • , C'2n-i'2n-i)^ ~ W(-^2n-i) denote a uniformly distributed 
vector on the moment space M.2n-i- From dsn) it follows that 

P {C^n-i e ^L-i) = 1 , 

and consequently the random canonical moments, say Pi 2n-i, • • • , -P2n-i,2n-i, corresponding to 
the random vector C2n-i are well defined with probability 1. Observing the representation (12.71) 
it follows that the density of the random vector (-Pi,2n-i, • • • , -P2n-i,2n-i) is given by 

2n-2 _rf ■\ 2n-2 

p-i) 11(^1^ n(p.(i-p.))^"-'-'. 

which means that {(Pj,2n-i)j=i^ | n G N} is a triangular array of rowwise independent random vari- 
ables, where P;,2n-i has a symmetric Beta-distribution on the interval [0, 1] with parameter 2n — j, 
that is P;,2n-i ~ B{2n — i, 2n — j). In what follows, let Pm,n{x) denote the mth random monic or- 
thogonal polynomial associated with the random vector C2„_i by equation (12.21) and Xi . . . , Xm,„ 
the corresponding roots which are real with probability 1 [see Szego (1975)]. Our first result gives 
an explicit representation for the joint density of the random vector (Xi_„, . . . , X„ „)-^. 

Theorem 3.1. The joint density of the roots Xi^n, • • • , -^n,n of the monic random orthogonal poly- 
nomial Pn,n{x) corresponding to a random vector C2n-i = (C'i,2n-i; ■ ■ ■ , C2n-i,2n~i)^ ~ W(-^2n-i) 
by equation 112. 2^) is given by 

n 

(3.2) f{Xi,...,Xn)=C- Yl Xj\'^Yiho^](^i)^ 

l<i<j<n i=l 

where the normalizing constant c is defined by 

1 r(2r + 2n) 

^ ~ r(n + 1) j-j^ r(2r + l)2r(2r + 2)' 
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Proof. Consider the random canonical moments (Pi 2n-i, • • • , -P2n-i,2n-i)^ corresponding to the 



random vector (C 



l,2n-l, 
'I -P. 



, C*2n-l,2n-l) 



2n-l 



Let 



\2n-1 = A,2n-i and for j >2 



i-l,2n- 



j,2n-l, 



then it follows from the recursive relation (12.6^ that the 
random orthogonal polynomial Pn,nix) can be represented as the determinant of a symmetric 
tridiagonal matrix, that is 



(3.3) 



X 



~V ^2n-3^2n-2 
~V ^2n-3^2n-2 X — ^2n-2 — ^2n-l 

Consequently, the roots Xi „, . . . , X„_„ of the polynomial P„_„(x) are the eigenvalues of the random 
Jacobi matrix Jn n where for m < n 



( ^ 



(3.4) 



J, 



m,n 



V ^1^2 
V ^1^2 ^2 + ^3 



'2m-3--2m-2 



V 



2m-3'='2m--2 '='2m-2 + ^2m-l j 



Now define for j = 0, . . . , 2n — 1 the quantities = 2Pj+i^2n-i — 1 then ao, . . . , a2n-i are indepen- 
dent random variables and aj follows a symmetric Beta-distribution on the interval [—1, 1] with 
parameter 2n—j — l. Moreover, a straightforward calculation shows that the matrix J„ = AJn^n—In 
can be represented as 



(3.5) 



with entries 



J n 



( h ai 

ai 62 

V 



bk+1 



(1 - a2k-l)a2k - (1 + tt2fc-l)«2fc-2 
{(1 - a2fe-l)(l - tt2fc)(l + "2fc+l)} 



1/2 



and a_i = —1. This matrix has been considered recently by Killip and Nenciu (2004) in a 
more general context. By the results of these authors we obtain that the density of the random 
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eigenvalues Ai, . . . , A„ of the matrix Jn is given by c]^^^-^^^^ |Aj — Aj |^, where c is an appropriate 
normalizing constant. Transferring this result to the eigenvalues of the matrix Jn,n it follows that 
the density of the roots of the random polynomial Pn,n{x) is given by (13.2^ . □ 



The density defined by (13.2^ is a special case of the Jacobi ensemble in the symplectic case [see 
e.g. Mehta (2004)] which has found considerable interest in the recent literature [see e.g. Collins 
(2005) or Johnstone (2008) among others]. The density of the general Jacobi /^-ensemble is given 

by 

n 

l<i<:j<:n i=l 

where a„,6„ > —1, c"^ is a normalizing constant and /? > [see Killip and Nenciu (2004)]. 
While most authors consider the case where n — > oo, a„/ri a; bn/n P, much less attention 
has been paid to the case where a„ and bn are fixed. We use the explicit representation for the 
density of the roots of P„ „(x) in terms of eigenvalues of the random matrix J„ „ defined in (13.41) 
to derive asymptotic properties for the empirical distribution function of the random variables 
Y Y 



Theorem 3.2. Let 



1 " 

Fn{x) = -Y,I{Xj,n<x} 



n 
i=i 

denote the empirical distribution function of the roots of the monic random orthogonal polynomial 
Pn,nix) associated with the random vector C2n-i — (C*i,2n-i; • • • ; C*2n-i,2n-i)"^ ~ 2n-i), then 
for all X G [0, 1] 

lim Fn{x) = - a.s. 
TT Jo ^t{l - t) 



Proof. Define 



(3.6) 



( i 



2 
1 

2\f2 



\ 



1 

2v^ 

1 

2 

1 

4 



then a straightforward calculation shows that the characteristic polynomial of the matrix D^^n is 
given by the monic Chebyshev polynomial of the first kind on the interval [0, 1], that is 



(3.7) 



T„(x) := det(a;/n - -D„,, 



)-2n+l 



cos(n arccos(22; — 1)), 
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which has roots 



cos((2A;-l)^) + l 



(3.8) Xk,n = ; fc = l,...,n. 

Moreover, define Tk{x) = 2"^^-^ Tfc(x) for > 1, To(x) = l/y/2 and T{x) = (To(a;), . . . ,T„_i(x))^, 
then it follows from the relation cos((ri+l)z)+cos((n— l)z) = 2 cos(n2;) cos(2;) by a straightforward 
calculation that 

(3.9) Xk,n T{Xk,n) = Dn,n T{Xk,n) ] k = l,...,n. 

This shows that the vectors 

/ 2k — 1 2k — 1 

(3.10) tk,n = T{xk,n) = ( I/V2, cos( tt), . . . , cos((n - 1) tt) j ; k = l,...,n 

are the eigenvectors of the matrix D„ „ corresponding to the eigenvalues xi^n, ■ ■ ■ , Xn,n, respectively. 
In the following discussion define 



1 " 

Nn{x) = - ^ I{Xj^n < X} 



n 
i=i 

as the empirical distribution function of the roots of T„(a;) (or T„(a;)), then an elementary calcu- 
lation shows that for all x G [0, 1] 

(3.11) lim N^{x) = - [ -j^=. 

From Bai (1999) we have for the Levy distance between the empirical distribution functions F„ 
and Nr, the estimate 



1 " 1 

L^{Fn, Nn) < - ^ I - f <-tr {Jn,n - Dn,nf 

^ j=l ^ 

(3.12) < - V(Jn- A«)' + - V(Jn+i- Am)', 

i=l 1=1 

where < ■ ■ ■ < -^{n),n and a;(i),„ < • ■ ■ < a;(n),n denote the ordered roots of the polynomials 

Pn,n{x) and r„(x), respectively, and Jij and Aj are the elements of the tridiagonal matrices 
Jn,n and defined in (13. 4p and (13.61) . respectively. Using the notation Qj^2n-i = 1 — Pj,2n-i 
(j = 1, . . . , 2n — 1) it now follows by a straightforward calculation that for 2 < i < n 

{Jii — Dii)^ = (^2i-3,2n-lP2i-2,2n~l + Q2i-2,2n~lP2i-l,2n-l ~ 

r, 1, 1, 1| 1^ 1| 1|D 1 

^ j |v2j-3,2n-l ~ 2' ~'~ 2I ~ 2 |'^2i-2,2n-l ~ 2 ' ~'~ 2I ~ 2 

I -^|2| 1|2| 1|2 

^ |-P2i-3,2n-l 2I + 9|-P2j-2,2n-l — 2 I + | Aj-l,2n-l " 2 I 
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and (Jii — Dii) < \Pi^2n-i — 2I • Similarly we obtain for 2 < i < n — 1 

I |2 I 1 I 

\Jii+l — -Dji+l < Q2i-2,2n-l-P2i-l,2n-lQ2i-l,2n-l-P2i,2n.-l ^ 77: 
III I 

~ g |2|P2i-2,2n-l — 2I + 3|-P2j-l,2n-l ~ 2 I ~'~ l'^2i,2n-l " -^l 

and |Ji,2 — Di^2f < ||-Pi,2n-i ~ H + i|-P2,2n-i ~ ^ I • In tliG following discussion we will show that 



n-l 

2 a.s. 



(3.13) ^J] 1^^2.2^-1-^ 

i=l 

(3.14) |P2.-1,2„-1 - ^ 



0, 



2 a., 



1=1 



which directly implies 



1 " 



2 



^ 0. 



77, 

i=l 



Moreover, for the remaining sum in (13.120 it follows that 

-5^(J..+1-Am)' < 7 -El^2.2n-1--| + -El^2.-l,2n-l--| 

i=l i=l i=l 

= o(l), 

where the constant 7 does not depend on n. The assertion is now a consequence of (13.111) and 
(I3.12p . which yields for the Levy-distance between the empirical distribution function F„ and 
distribution function F of the arcsine measure 

L(F„,F)<L(F„,iV„) + L(iV,,F) = o(l) a.s. 

For a proof of the almost sure convergence in (13.131) and (I3.14p we restrict ourselves to the statement 
(I3.13P , the remaining case is treated similarly. In order to prove (|3.13p we will use a strong law of 
large numbers for arrays of rowwise independent random variables. To be precise, define 

Zn,k = (P2.,2„-i - 1? - g^2^ _ ^ 4 k = l,...,n-l 

then E[Zn^k] = and a tedious calculation shows that for A; = 1, . . . , n — 1 and any t > 

P{\Zn,k\ >t)< P{\Zn,n-l\ > t) , 

where the distribution of the random variable Zn,n-i does not depend on n [note that P2n-2,2n-i ~ 
5(2,2)]. Consequently, we obtain from Theorem 2 in Hu, Moricz and Taylor (1989) that 

n-l 



fc=i 



a.s. 



for any p G [1,2). Observing that 



n-l 



1 1 

4E 2(2n-l-2A:) + l =^(^"g"^ 



it follows for any p E [1,2) that 



ii:iA«.-.4r-;^^+«(^)}-( 

I 1 N 



n 



i-i/p 



a.s. 



k=l 

which establishes (13.131) and completes the proof of Theorem 3.2. □ 

Our final result refers to the asymptotic behaviour of the roots of the mth orthogonal polynomial 
associated with the random moment vector 6*2^-1 = (C*i,2n-i, ■ • ■ , C*2n-i,2n-i)^ ~ W(-^2n-i), where 
m is fixed. In this case, the limiting distribution is normal, where the roots of the Chebyshev 
polynomial T„(x) defined in (13. 7p are used for the centering. 

Theorem 3.3. Let Xi^n, ■ ■ ■ ,Xm,n denote the roots of the m-th random monic orthogonal poly- 
nomial Pm.,n{x), which corresponds to the random vector C2n-i ~ W(-^2n-i) by equation 112. 
then 

where xi^m, ■ ■ ■ ,Xm,m are the roots of the Chebyshev polynomial of the first kind defined by ( fg. 8\) 
and the matrix T is given by T = ■^{lk,i)'k'i=iJ where 

^ ^ m ^ jn— 1 

lk,i = 4 + 2 XI '^j~i(^k,m) Tf_^{xi,m) - 4 5^ (r/_i(xfe,™) T/(a;i,™) + T/_i(x«,^) T/(xfc,„)) 

J=2 j=l 



+ ^ Ti{Xk,m) Ti{xi^rn) + " Tj^i{Xk,m) Tj{Xk^m) Tj_i{xi^rn) Tj{xi^rn) 

i=2 

^ m—2 

~ 7 ^ ] (Tj^i(Xk^m) Tj(Xk^m) Tj^Xi^m) Tj^i(xi^m) ~^ Tj_i(xi^m) Tj^Xi^m) Tj(Xk,m) Tj^i(Xk,m)) 



4 

i=l 



Tj{x) = cos(j arccos(2a; — 1)) and Xk,m is the kth zero of the polynomial Tm{x) defined in 113. 8\) . 
Proof. It is easy to see that for fixed m G N we have for = 1, . . . , m 
(3.15) {PkM^i - I) ^AfiO,l), 
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which implies that the vector 4:^/n{Pl 2n-l, ■ ■ ■ , Pm,2n-i)'^ is asymptotically multivariate normal 
distributed with mean and covariance matrix Im (note that Pi 2n-i, . . . , Pm,2n-i are independent 
random variables). A straightforward calculation now shows that 



(3.16) Ay/n I ^Si,2n-l, • • • , S2m-l,2n-l 

where the matrix A^m-i G M2m-ix2m-i -g tridiagonal and given by 



2' 4'' "'4/ J '^=^ 



A/'2„_i(0,A 



2m-\j 1 



( 1 



(3.17) 



1 
'2 

1 _1 

2 4 

1 1 _1 

"4 2 4 



V 

A further application of the delta-method finally yields 



_ 1 

4 

1 1 



(3.18) 



4yri 



^2 -r --3 



-2m-2 + ^2m-\ 



'\^2 



2m-3^2m-2 



\ f 1/2 \ 

1/2 

1/2 
l/(2v^) 
1/2 

J \ J 



V 



where the covariance matrix V = diag{Vu, V22) € M^"* ix2m 1 jg block diagonal with blocks 
= v4„ e M'"^'" and 



/ 1 -7S 

4. 1 



22 



\ 



V2 



.1 1 _1 

2 2 



V 



nm— Ixm— 1 



1 / 



Consequently, the difference of the matrices Jm,n — Dm,n defined by (13.41) and (13.61) converges also 
weakly, that is 

lim 4v^ {Jm,n - Dm) = S (weakly). 
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where S denotes a triangular matrix given by 



S 



Ni Ma N2 

N2 M3 A^3 



V 



Nm-l Mm j 



and {Ml, . . . , M^} and {A^i, . . . , N^-x} define independent samples Mi ~ A/'(0, 1), M, ~ A/'(0, 1/2) 
(j = 2, . . . , m) ~ A/'(0, 1/8) (j = 1, . . . , m - 1) such that 



(3.19) 



(3.20) 



Gov (Mi,Mj) 



Gov (Ni,Nj) 



— \ if z = l,j = 2 or 2 = 2,j = 1 
-i if j = ^ + 1 or j = i - 1 
else 

if z = l,J=2 or z = 2,j = l 
if j = i + 1 or j = i - 1 
else 



Using the same arguments as in Dimitriu and Edelman (2005) it now follows that the first order 
properties of the matrix Jm,n are the same as the first order properties of the matrix -D^.n + 4^*^- 
In particular, using Lemma 2.1 in this reference we obtain for i = 1, . . . ,m that the asymptotic 
distribution of the vector 



coincides with the distribution of the random vector 



^,^^) } 



(3.21) 



G 



J.T C i- y^T C J- 



t 



where the vectors tim, ■ ■ ■ ,t 



, . . . , um,in 



are defined by (13.10p . Obviously, G has a multivariate 
normal distribution with mean zero and it remains to calculate the corresponding covariance 
matrix. For this purpose we note that for k = 1, . . . , m 



m— 1 



2A; — 1 m 
2m ' 2 



Therefore it remains to calculate the covariance matrix of the vector (tf ^ ^ ti,m, ■ ■ ■ ^tmm tm,m 
For this purpose we denote by (t>i, . . . , Vm) and {wi, . . . , Wm) the components of the vectors tk- 
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and ti^m, respectively, and obtain 



m m m—1 



m—1 



Gov {tl^ S h,m, S ti,m) = Coy(^J2''^^^^J2'^J^^) +4 Gov 



viVi+iNi, 

i=i 



i=l j=l i=l 



m _ m—1 




+ <i^')} 




m-2 



1=1 



ViVi+iWiWij^l - - 



V2W2W3 + W2V2V3 + > ( 



ViVi+iWi+iWi+2 + Vi+iVi+2WiWi+i 



i=2 



where we have used (13.190 and (13.201) and wi = vi = l/\/2. This proves the assertion of the 
Theorem. □. 
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